Decoherent Histories Analysis of Minisuperspace 
Quantum Cosmology 

J.J.Halliwell 

Blackett Laboratory 
Imperial College 
London SW7 2BZ 
UK 

Abstract. 

Recent results on the decoherent histories quantization of simple cosmological models 
(minisuperspace models) are described. The most important issue is the construction, from 
the wave function, of a probability distribution answering various questions of physical interest, 
such as the probability of the system entering a given region of configuration space at any stage 
in its entire history. A standard but heuristic procedure is to use the flux of (components of) 
the wave function in a WKB approximation as the probability. This gives sensible semiclassical 
results but lacks an underlying operator formalism. Here, we supply the underlying formalism by 
deriving probability distributions linked to the Wheeler-DeWitt equation using the decoherent 
histories approach to quantum theory, building on the generalized quantum mechanics formalism 
developed by Hartle. The key step is the construction of class operators characterizing questions 
of physical interest. Taking advantage of a recent decoherent histories analysis of the arrival time 
problem in non-relativistic quantum mechanics, we show that the appropriate class operators 
in quantum cosmology are readily constructed using a complex potential. The class operator 
for not entering a region of configuration space is given by the S-matrix for scattering off a 
complex potential localized in that region. We thus derive the class operators for entering one 
or more regions in configuration space. The class operators commute with the Ifamiltonian, 
have a sensible classical limit and are closely related to an intersection number operator. The 
corresponding probabilities coincide, in a semiclassical approximation, with standard heuristic 
procedures. 



1. Introduction 

1.1. Opening Remarks 

The problem of finding a sensible quantization of the Wheeler-DeWitt equation of 
minisuperspace quantum cosmology, 

H-^ = Q (1) 

continues to attract considerable interest [1, 2, 3, 4, 5]. Although the setting of this problem 
is simple cosmological models with just a handful of homogeneous parameters, the techniques 
employed in answering this question may be relevant to general approaches to quantum gravity, 
such as the loop variables approach or causal set approach. This is because the central difficulty 
in consistently quantizing and interpreting the Wheeler-DeWitt equation is the absence of a 
variable to play the role of time and all approaches to quantum gravity must confront this issue 
at some stage [6]. 



A frequently studied example consists of a closed FRW cosmology with scale factor a = e° 
and a homogeneous scalar field cf) with (inflationary) potential ^(0) [7]. The Wheeler-DeWitt 
equation for this model is 



Given suitable boundary conditions, one can solve this equation for the wave function ^'(a, 0) 
and attempt to use it to answer a number of interesting cosmological questions. There are many 
such questions: Is there a regime in which the universe behaves approximately classically? What 
is the probability that the universe expands beyond a given size ao? What is the probability 
that the universe has a certain energy density at a given value of the scale factor? What is the 
probability that the universe's history passes through a given region A of configuration space, 
characterized by certain ranges of a and 0? 

Such questions, of necessity, do not involve the specification of an external time. Classically, 
the system's trajectories in minisuperspace may be written as paths {a{t) , (j){t)) , but here t is 
a convenient but unphysical parameter that labels the points along the paths - it does not 
correspond to the physical time measured by an external clock. The absence of a physical 
time is reflected in the quantum theory by the fact that the quantum state obeys the Wheeler- 
DcWitt equation Eq.(l), not a Schrodingcr equation, and it is this difference that presents such 
a challenge to conventional quantization methods. 

Although very plausible heuristic semiclassical methods exist for formulating and answering 
the above questions (in particular, the WKB interpretation [1]), it is of interest to see whether 
or not there exists a precise and well-defined quantum-mechanical scheme underlying these 
heuristic methods. We are not looking for high standards of mathematical rigour - just a 
standard quantum-mechanical framework of operators, inner product structures etc, obeying 
some reasonable requirements. A broad framework along these lines is the generalized quantum 
mechanics of Hartle [8] . The purpose of this paper is to carry out a specific implementation of 
Hartle's general framework, building also on earlier related attempts [9, 10, 11, 12, 13, 14, 15]. 
A much more detailed account of this work is given in Ref . [16] . 

1.2. Inner Products and Operators for the Wheeler-DeWitt Equation 

In the Wheeler-DeWitt equation (1), H is the Hamiltonian operator of a minisuperspace model 
with n coordinates g", and is typically of the form 



where V is the Laplacian in the minisuperspace metric fab, which has signature ( — !- + + •••)• 
It is naturally linked to the current. 




H=-V'^ + U{q) 



(3) 




(4) 



which is conserved 



V- J = 



(5) 



Closely associated is the Klein-Gordon inner product defined on a surface S 




(6) 



where da"' is a surface element normal to S. In flat space with a constant potential, the Wheeler- 
DeWitt equation is just the Klein-Gordon equation. Its solutions may be sorted out into positive 



and negative frequency in the usual way. With a httle attention to sign, it is then possible to 
use components of the current to define probabilities. 

However, in general, it is not possible to sort the solutions to the Wheelcr-DeWitt equation 
into positive and negative frequency. This is one manifestation of the problem of time and more 
elaborate methods are required to associate probabilities with the Wheeler-DeWitt equation. 
There arc two main issues: finding an inner product, and then finding suitable operators. 

The issue of finding a suitable positive inner product is reasonably straightforward and goes 
by the name of Rieffel induction, or the induced (or physical) inner product [9, 17]. We consider 
first the usual Schrodinger inner product, 

(*i|^2)= J dy^l{q)^2{q) (7) 
We then consider eigenvalues of the Wheeler-DeWitt operator 

H\^xk) = AI^Afc) (8) 

where k labels the degeneracy. These cigcnstatcs will satisfy 

{^x'k'\'^xk)=S{\-X')6ik-k') (9) 

from which it is clear that this inner product diverges when A = A'. The induced inner product on 
a set of eigenstates of fixed A is defined, loosely speaking, by discarding the (5-function 6{X — A'). 
That is, the induced (or physical) inner product is then defined by 

{^Xk'\^Xk)phys = Sik-k') (10) 

This procedure can be defined quite rigorously, and has been discussed at some length in 
Refs.[17, 9]. It is readily shown that the induced inner product coincides with the Klein-Gordon 
inner product when a division into positive and negative frequencies is possible, with the signs 
adjusted to make it positive. (This is described in Appendix A of Ref.[16]). 

Turning now to the construction of interesting operators, the interesting dynamical variables 
associated with the Wheeler-DeWitt equation are those that commute with H. This is because 
the constraint equation is related to reparametrization invariance - which is reflected in the 
absence of a physical time variable - and we seek operators which are invariant. A wide class of 
operators commuting with H are of the form 

/oo 
dt B{t) (11) 
-oo 

which clearly commute with H since 

/oo 
dt B{t + s) 
-oo 

= A (12) 

Many examples are given in Refs.[5, 17, 11]. However, another way of constructing such operators 
involves taking products, 

oo 

^= n ^(*) (13) 

t = — QO 

which may be shown to commute with H using essentially the same argument [12], but clearly 
further mathematical detail is required to give meaning to the infinite product. (Here t is the 
unphysical parameter time). 



91 



Figure 1. Given a solution 'I' to the Wheeler-DeWitt equation, what is the probabihty that 
the system enters a series of regions Ai, A2 • • • in configuration space, or crosses a surface S, at 
any stage in the system's entire history? 

Given these prescriptions for inner products and operators, one may then attempt to construct 
operators and probabihties implementing some of the questions mentioned above. We will focus 
on the following general question: Given a solution ^ to the Wheeler-DeWitt equation, what 
is the probability of finding the system in a region A of configuration space, or of crossing a 
surface E, at any stage in the system's history? The question is depicted in Figure 1. 

Questions involving surface crossings are not unlike more familiar questions in non-relativistic 
quantum mechanics, where there is a physical time parameter, but the key difference in quantum 
cosmology is that even classically, a given trajectory will typically cross a fixed surface more 
than once. It is precisely these sorts of issues that need to be carefully phrased in the quantum 
theory. Whilst the operator formalism briefly outlined above has been used to address such 
questions [5], the problems of characterizing properties of trajectories and surface crossings in 
quantum theory is naturally accommodated in the decoherent histories approach to quantum 
theory. 

1.3. The Decoherent Histories Approach 

A more general approach to quantizing and interpreting the Wheeler-DeWitt equation is 
provided by the decoherent histories approach to quantum theory, and this approach will be 
the focus of this paper [18, 19, 20, 21, 22, 23, 24, 25, 26]. In this approach, probabilities are 
assigned to histories using the formula 



Here, p is the initial state (in our case a pure state) and Ca is a class operator characterizing 
the histories a of interest. In non-relativistic quantum mechanics these class operators are given 
by time-ordered strings of projection operators. 




(14) 



Ca = Pan (*n) ' ' ' Pa2 {t2)Pai (*l) 

(or by sums of such strings). The class operators always satisfies the condition 



(15) 




(16) 



a 



For the reparametrization invariant systems considered here, the definition of the class operators 
is more subtle, and we return to this below. 

Because of interference between pairs of histories, probabilities cannot always be assigned. 
To check for interference we therefore consider the decoherence functional, 

D{a,a') = Tr(^CapCl) (17) 

When 

D{a,a') = (18) 

for all pairs of histories in the set with a ^ a' , wc say that there is decoherence of the set of 
histories and probabilities may be assigned Eq.(14). When there is decoherence, it is easily seen 
from Eq.(16) that the probabilities Eq.(14) arc also given by 

p{a) = Tr{CaP) = Tr(cip) (19) 

Decoherence guarantees that this expression is real and positive, even though the class operators 
are not positive or hermitian operators in general. 

The structure of the decoherent histories approach is very general and may be applied to a 
wide variety of situations, given an initial state, class operators, and a suitable inner product 
structure with which to construct the decoherence functional. A useful formulation of the 
decoherent histories approach encapsulating this generality is Hartle's generalized quantum 
mechanics [8]. This in essence defines a class of quantum theories through a decoherence 
functional obeying some simple requirements but does not rely on the specific form Eq.(15) of the 
class operators used in non-rclativistic quantum mechanics. That framework is the background 
of what we do here, but it differs from the specific implementation of the framework in Ref.[8] 
(and elsewhere) through our choice of inner products and class operators. 

For the application to the Whcclcr-DeWitt system considered here, the initial state is take 
to be a solution to the Wheeler-DeWitt equation and the inner product is the induced inner 
product described above. The most crucial element is the specification of the class operators 
Cq,. As indicated, in generalized quantum mechanics they can be more general than the non- 
relativistic version, Eq.(15). The product in the string can be taken to be continuous time [26]. 
Also, it is often valuable, sometimes essential, to allow the projectors to be replaced by more 
general operators, such as POVMs. The class operators must also properly characterize the 
histories one is interested in. It is not always obvious how to do this but useful clues often come 
from looking at the classical analogue of the class operator (where all the projectors commute). 

Here, we are interested in histories which enter a region A of configuration space, or which 
cross a surface S, but without regard to time. This absence of a physical time variable seems 
particularly challenging, given that time seems to be central to the definition of non-relativistic 
analogue, Eq.(15). Closely related to this is the role of the constraint equation, Eq.(l). As 
noted already, these two features are directly related to the underlying symmetry of the theory 
- reparametrization invariance - and this symmetry is respected if the class operators commute 
with the constraint, 

[^,CJ=0 (20) 

Eq.(20) is in keeping with standard procedures of Dirac quantization [5]) and also ensures that 
the class operators have a sensible classical limit [12]. 

Some partially successful attempts to construct such class operators have been given 
previously [9, 10, 11, 12], but ran into various problems to do with the Zeno effect and 
with compatibility with the constraint equation. The main aim of this paper is to give fully 
satisfactory definitions of class operators for quantum cosmological models and explore their 
decoherence properties and probabilities. 



1.4- Some Properties of the Wheeler- De Witt Equation and the WKB Interpretation 
To prepare the way for the full decoherent histories analysis of quantum cosmology, it is 
important to discuss some properties of the Wheeler-DeWitt equation and review the commonly 
used heuristic semiclassical interpretation of the wave function, since a proper quantization must 
recover this structure in some limit. 

The Wheeler-DeWitt equation Eq.(3) is a Klein-Gordon equation in a curved configuration 
space with indefinite metric fabio) and potential U{q) which can be positive or negative. The 
curvature effects of the metric are not significant in relation to the issues addressed in this paper, 
so we will assume for simplicity that the metric is flat. 

The classical constraint equation corresponding to the Wheeler-DeWitt equation Eq.(3) is 

-JabQ^q'' + [/ = (21) 

from which one can see that the classical trajectories are timelike in the region J7 > and 
spacelike in ?7 < 0. (The timelike direction is that of increasing scale factor and the spacelike 
directions correspond to matter and anisotropic modes). The quantum case has analogous 
features. In simple models such as Eq.(2), the character of the solutions to the Wheeler-DeWitt 
equation depends on the sign of U . For large scale factors, C/ > and the wave function is 
oscillatory, corresponding, very loosely, to a quasiclassical regime, and for small scale factors, 
?7 < 0, and the wave function is exponential, corresponding to a classical forbidden regime. 
However, there are certain types of models (such as those with an exponential potential for the 
scalar field), in which the identification of the oscillatory and exponential regions depends also 
on whether the constant U surfaces are spacelike or timelikc [27]. We will not address this here. 
One can also associate a Feynman propagator Gp with the Wheeler-DeWitt operator, 

oo 

iHt-et 



-L 



Gf = / die 



—I 



H-ie 



(22) 



where e 0+. (Numerous propagator-like objects of this type have been considered in quantum 
cosmology [28]). Locally, on scales smaller than the scale on which the potential U significantly 
varies, the propagator Gp{q,q') will be essentially identical in its properties with the analogous 
object for the relativistic particle for flat space. Therefore, in the region [/ > 0, for points q 
and q' which are timelike separated, it will propagate positive frequency solutions to the future 
and negative frequency solutions to the past. It will be exponentially suppressed for initial and 
final points q and q' which are spacelike separated (this is the familiar "propagation outside the 
lightcone" effect). For ?7 < 0, similar statements hold but with timelike and spacelike reversed. 
Similar features hold on larger scales in a semiclassical approximation. These properties are 
important to understand the class operator constructed below. 

Very plausible but heuristic answers to questions concerning crossing surfaces and entering 
regions are readily found using the WKB approximate solutions to the Wheeler-DeWitt equation 
and the Klein-Gordon current [1]. In the oscillatory regime, the Wheeler-DeWitt equation may 
be solved using the WKB ansatz 

* = Re'^ (23) 
where the rapidly varying phase S obeys the Hamilton-Jacobi equation 

{VSf + J7 = (24) 

and the slowly varying prefactor R obeys 



V • (|i?|^V5) = 



(25) 



Figure 2. A WKB wave function with phase S corresponds to a set of classical trajectories 
with tangent vector VS. The probability for entering a region A is the amount flux of the wave 
function intersecting A. It may be expressed either in terms of the ingoing flux across Ej„ or 
equally, in terms of the outgoing flux across ^out- 



The latter equation is just current conservation for the WKB current, 

J = \R\'^VS (26) 

Wave functions of the WKB form Eq. (23) indicate a correlation between position and momenta 
of the form 

p = VS (27) 

and this suggests that the wave function Eq.(23) corresponds to an ensemble of classical 
trajectories satisfying Eq.(27). The current J may then be used to define a measure on this set 
of trajectories. For example, we consider a surface S and choose the normal n" to the surface 
so that n • VS > 0. Then the probability of the system crossing the surface is taken to be 

pi^)= J^d^^-^qn'^Ja (28) 

Of particular interest here is the probability of entering a region A. This is clearly related to 
the flux at the boundary of the region. The current will typically intersect the boundary of A 
twice. However, we can split the boundary S of A into two sections: Sj^ at which the current 
is ingoing and Hout at which the current is outgoing. The probability of entering A is may then 
be expressed in the two forms 

P(A) = - / d"-^(?nV„ 
= / rf-Ja 

where here we have defined the normal n" to point outwards. See Figure 2. These two forms 
are equivalent since the current is locally conserved. 

Typically, little is said of the regions in which the wave function is exponential, except that 
they are similar to tunneling regions in non-relativistic quantum mechanics, so are classically 



(29) 



forbidden in some sense. However, the wave function is not necessarily small in these regions, so 
there is surely more to it than this. In Ref.[29], it was noted that, unlike the oscillatory regions, 
the exponential regions do not indicate a correlation between position and momenta of the form 
Eq.(27), and it seems that this should be significant in some way. 

Whilst the WKB interpretation is very plausible and adequate for most situations of interest, 
it leaves many questions unanswered. The main issue is to understand the operator origin of 
these probabilities, in terms of the language of operators commuting with H developed above. 
Furthermore, what can one say about superpositions of WKB states? Can interference between 
them be neglected? Also, what can one say about the exponential, classically forbidden regions? 
Is there a more precise way of saying that they are non-classical? 

1.5. This Paper 

The purpose of this paper is to present a decoherent histories quantization of the Wheeler- 
DeWitt equation, and in particular to exhibit exactly defined class operators which characterize 
histories entering regions of configuration space or crossing surfaces, without reference to an 
external time. The key idea is to use a complex potential to define the class operator for not 
entering a region A in configuration space. In particular, we take the class operator for not 
entering to be the S matrix describing scattering off a complex potential localized in A. This 
turns out to have all the right properties - it avoids the Zeno effect, is compatible with the 
constraint and has a sensible classical limit. 

In Section 2, to explain and motivate the use of a complex potential we review the use 
of such potentials in the decoherent histories analysis of the arrival time problem in non- 
relativistic quantum mechanics. In Section 3, we describe the construction of class operators for 
the Wheelcr-DcWitt equation using a complex potential. The properties of the class operator 
for entering a region A are described in Section 4. The class operator has a sensible classical 
limit which, crucially, registers just one intersection of a trajectory with a surface, even when 
the trajectory intersects twice or more. In the quantum case, the class operator is an operator 
describing the ingoing flux across the boundary of the region, an expected result on semiclassical 
grounds, and is closely related to an intersection number operator. 

In Section 5 we consider the WKB regime and show that the decoherent histories analysis 
reproduces the expected heuristic interpretation of the wave function. We summarize and 
conclude in Section 6. 

2. The Arrival Time Problem in Non-Relativistic Quantum Theory 

In this section we summarize some of the key features of the arrival time problem in non- 
relativistic qiiantum mechanics [30]. These details arc very relevant to the quantum cosmology 
case and in particular motivate the use of complex potentials in the definition of class operators. 

In the one-dimensional statement of the arrival time problem, one considers an initial wave 
function \ip) concentrated in the region x > and consisting entirely of negative momenta. The 
question is then to find the probability n(r)(ir that the particle crosses x = between time r 
and r -|- dr. See Figure 3. 

The canonical answer is the current density 

J{r) = ^(V'rl + 6{x)p) IV'r) (30) 

where I'i/'-r) is the freely evolved state. This has the correct semiclassical limit, but can be 
negative for certain types of states consisting of superpositions of different momenta (backflow 
states). It is of interest to explore whether this simple semiclassical result emerges from more 
elaborate measurement-based or axiomatic schemes. Many such schemes naturally lead to the 
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Figure 3. The quantum arrival time problem in non-relativistic quantum mechanics. Given an 
initial state localized entirely in x > and consisting entirely of negative momenta, what is the 
probability that the particle crosses the origin during the time interval [ti,t2]? 

use of a complex potential —iV{x) in the Schrodinger equation [31], where 



where H is the free Hamiltonian. The idea here is that the part of the wave packet that reaches 
the origin during the time interval [0, r] should be absorbed, so that ('0(r)|^/'(r)) is the probability 
of not crossing x = during the time interval [0,t]. The probability of crossing between r and 
r + dr is then 



This may be approximately evaluated in the limit Vq <^ E (where E is the typical energy scale) , 
which is the limit of negligible reflection off the complex potential. The result is 

n(r) = 2Vo r dt e-2^o(^-*) J{t) (34) 
The probability for crossing during a finite interval [ri , T2\ is then given by 



v{x) = v^ei-x) 



(31) 



With such a potential, the state at time r is 



IV'(r)) = exp {-iHr - Vb^(-a;)r) 



(32) 



n(r) = --(V'(r)|V'(r)) 



(33) 




n 



(35) 



and if this time interval is sufficiently large compared to I/Vq, we have 




Tl 



(36) 



which means that the dependence on the potential drops out entirely at sufficiently coarse grained 
scales, and there is agreement with semiclassical expectations involving the current J{t). 

In the decoherent histories analysis of the arrival time problem [32, 33], we first consider the 
construction of the class operator Cnc for not crossing the origin during the finite time interval 
[0,r]. We split the time interval into N parts of size e, and the class operator is provisionally 
defined by 

Cnc = Pe-'^'P ■ ■ ■ e-'^'P (37) 

where there are N + 1 projections P = 9{x) onto the positive x-axis and N unitary evolution 
operators in between. One might be tempted to take the limit AT — > oo and e ^ 0, but this 
yields physically unreasonable results. This limit actually yields the restricted propagator in 
x>0, 

Cnc = gr{T,0) 

= Pe-'^"^^ (38) 
This object is also given by the path integral expression 

{xi\gr{T,0)\xo) = J^Vxexp{iS) (39) 

where the integral is over all paths from x{0) = xq to x(r) = xi that always remain in x{t) > 0. 
However, the class operator defined by Eq.(38) has a problem with the Zeno effect - it consists 
of continual projections onto the region a; > and as a result the wave function never leaves the 
region. This is refiected in the fact that the restricted propagator is unitary in the Hilbert 
space of states with support only in a; > 0. This is a serious difficulty which has plagued a 
number of earlier works in this area [34, 35, 36]. 

To avoid the Zeno effect, the key is to keep e non-zero. The class operator Eq.(37) is not 
easy to work with for finite e, but fortunately a result first suggested by Echanobe et al. comes 
to the rescue [37, 38]. This is that the string of operators in Eq.(37) is in fact approximately 
equivalent to evolution in the presence of the complex potential —iV introduced above. That 
is, 

Pe-'^'P---e-'^'P « ei^p{-iHoT -Vo9{-x)t) (40) 

This connection is valid as long as e <C 1/{AH) and e and Vq are related by Vbe fa 4/3 [38]. 
It strongly suggests that the class operator Cnc, normally defined by a string of projection 
operators, is justifiably defined instead using a complex potential. That is, we define 

Cnc = exp i-iHoT - Voei-x)T) (41) 

The subsequent decoherent histories analysis was described in detail in Ref.[32]. The 
corresponding class operator for crossing during a time interval [ri , T2] was shown to be 

Cc{t2,ti)= r' e-^^°("-*Ve-*-^°*-^* (42) 

and again in the approximation Vq <^ E and for time intervals greater than I/Vq this may be 
shown to have the simple and appealing form 



Tl 



T2 „ 

2m 



Cc{r2,n)^ I dte-^''-'-^ {p6{xt) + d{xt)p) (43) 



which is now independent of the Vq. (These definitions of class operators differ by a factor 
of exp(— zi^r) from those defined in Section IC, but this does not make any difference in the 



decoherence functional and probabilities). With this class operator, one can show that there is 
decoherence of histories for a variety of interesting initial states, such as wave packets (but not 
for superposition states with backflow), and for such states, the general result Eq.(19) implies 
that the probability for crossing is simply 

p(r2,Ti) = (*|Cc(r2,Ti)|*) (44) 

which agrees precisely with the semiclassically expected result, Eq.(36). 

In summary, the dccohercnt histories analysis of the arrival time problem in non-relativistic 
quantum mechanics indicates that it is reasonable to define class operators for not entering a 
space time region using a complex potential, as in Eq.(41), and that such a definition gives 
sensible a semiclassical limit, independent of the potential, at sufficiently coarse grained scales. 

3. Construction of the Class Operators Using a Complex Potential 

We now come to the central issue concerning this paper, which is the construction of class 
operators for the decoherent histories analysis of the Wheeler-DeWitt equation. We seek class 
operators for the system Eq.(3) describing histories which enter or do not enter the region A, 
without specification of the time at which they enter. It is easiest to first focus on the class 
operator Ca for not entering and the class operator for entering is then given by 

Ca = 1 - Ca (45) 

The earliest attempts to define class operators for the Wheeler-DeWitt equation involved 
defining C\ as a sort of propagator obtained by integrating restricted propagators of the form 
Eq.(39) over an infinite range t (now regarded as the unphysical parameter time) [9, 10, 11]. 
However, in addition to having problems with the Zeno effect (which was not in fact appreciated 
in these earlier works), such constructions are difficult to reconcile with the constraint equation 
and some ad hoc modifications of the basic construction were required to give sensible answers. 

A rather different approach to constructing Ca was given in Ref.[12]. This was again 
problematic, but we review the construction here, since it is readily modified to yield a STiccessful 
definition of the class operators. Wc denote by P the projector onto A and P the projector onto 
the outside of A. Our provisional proposal for the class operator for trajectories not entering A 
is the time-ordered infinite product, 

oo 

= n ^(*) (46) 

t=— oo 

where t is the unphysical parameter time. Subject to a more precise definition, given shortly, 
this object has the required properties. It is a string of projectors. Classically, it is equal to one 
for trajectories which remain outside A at every moment of parameter time. Also, it commutes 
with H, at least formally. 

To define this more precisely, we first consider the product of projectors at a discrete set of 
times, ii, ^1 + e, ii + 2e, • • • , ii + ne = t2- We define the intermediate quantity, C/\{t2, h) as the 
continuum limit of the product of projectors, 

CA{t2, h) = lim P{t2)P{t2 - e) . . . P{ti + e)P{ti) (47) 

e— ^0 

where the limit is n ^ oo, e — )■ with t2 — h fixed. The desired class operator is then 



Ca= lim CA{t2,ti) 

t2— >-00,tl— > — oo 



(48) 



The class operator operator is clearly closely related to the restricted propagator (the 
generalization of Eq.(39)) in the region outside A, since we have 



CAit2,ti) = e^^*2 grit2,ti) e-^^*i (49) 

and therefore 

(7a = lim e'"^^ gr{t2M) (50) 
The class operator commutes with H. This is because, from Eq.(49) 

e'^'CA{t2,ti)e~'^" = e'^^^^+'^ gr{t2, h) e-^^(*i+^) (51) 
This becomes independent of s as t2 — >■ oo, ti —oo, hence 

[H,Ca] = (52) 

The problem with this definition, however, is that it suffers from the Zeno effect, exactly like 
the analogous expression Eq. (37) (in the limit e — >■ 0) in the non-relativistic arrival time problem. 
But the key idea here is that we may also get around the problem in the same way, using a 
complex potential. That is, we "soften" the restricted propagator and make the replacement 

9r{t2,h) ^ exp {-iH{t2 - h) - V{t2 - h)) (53) 

where V{q) = Vo/aI?)- Here, Vq > is a constant and /a (9) is the characteristic function of 
A, so is 1 in A and outside it. Or we may equivalently write V = VqP, where recall, P is the 
projector onto A. This means that our new definition for the class operator is 

Ca= lim e'"^^ exp{-i{H -iV){t2-ti)) e-'^^' (54) 

The class operator thus defined has an appealing form: it is the S-matrix for system with 
Hamiltonian Eq.(3) scattering off the complex potential V. As required, it commutes with H. 
This is the most important definition of the paper. 

Now an important observation. The class operator derived above has been defined as a 
time-ordered product of operators in which the direction of parameter time increases from right 
to left. However, since parameter time is unphysical, there is absolutely no reason why the 
parametrization should not run in the opposite direction. This produces an operator which is 
the hermitian conjugate of Eq.(54). As one can see from Eq.(19), this makes no difference in 
the final expressions for probabilities. In fact, it is most natural to define the class operator in 
such a way that it is invariant under reversing the direction of parametrization. We thus define 
a modified class operator which is hermitian: 

C'A = \{CA + Ci) (55) 

In what follows, for simplicity, we will primarily work with the non- hermitian class operator 
Eq.(54) and revert to the hermitian one, Eq.(55) where appropriate. The difference between 
them will turn out to be significant only for the class operator for two or more regions. 

We now cast the above class operator in a more useable form. The following identities are 
readily derived: 

^-i{H-iV){t2-ti) ^ ^-iH{t2-ti) _ ^-iH{t2-t)y^-i{H-iV){t-tl) ^5g^ 

Jtl 

r-t2 



H(t2-tl) _ ^-i{H-iV){t2-t)y^-iH{t-ti) (^57^ 



Inserting the second expression in the first, we obtain 

-iiH-iV)it2-ti) ^ ^-iH(t2-ti) _ f^"" ^-iH{t2-t)y^-iH{t-tl) 
'■t2 ft 



tl 



+ r dt f ds e~'H{t2~t)y^-i{H-iV){t-s)y^-iH{s-tl) (58^ 



Inserting in the expression for the class operator (54) and taking the hmit, we obtain 



/oo 
dt V{t) 
-OO 

/OO pt 
dt / ds e^"ty^-iiH-iV){t-s)y^-iHs (59) 
-oo J — oo 

The class operator for entering the region is therefore given by 

/oo roo rt 

dt V{t) - dt ds eiHty^-i{H-iV){t-s)y^-iHs (gQ) 

-oo J —oo J — oo 

This is an exact and useful form for the class operator for entering A, and it is easily confirmed 
that it is of the form Eq.(ll) so commutes with H. 

Now we use a simple but useful semiclassical approximation. Noting that V = VqP, where 
P recall projects into A, note that the expression 

y^-i{H-iV){t-s)y (gj) 

describes propagation with the complex Hamiltonian H — iV between two points that lie inside 
A. This can easily be represented by a path integral in which it seems plausible that the 
dominant paths between this two end-points will lie entirely inside A, as long as the boundary 
is reasonable smooth and the semiclassical paths are not too irregular. If this is true, we may 
replace V = VoP by the constant complex potential V = Vo, that is, 

y^-i{H-iV)(t-s)y _ y^-iiH-iVo)it-s)y (•g2) 

Propagation with a complex potential in Eq.(61) will also involve reflection off the boundary of 
the region, in which case the semiclassical approximation Eq.(62) may fail. However, reflection 
is small for sufficiently small Vq [32, 31], and we will see this in more detail in Section 5. Hence 
we expect the semiclassical approximation Eq.(62) to hold for small Vq. 

With this useful approximation (which does not affect the fact that the class operator 
commutes with H), we have 

/oo poo rt 

dt V{t) - dt ds V{t)V{s)e-^°^*-''^ (63) 
-oo J—oo J —oo 

Again using V = VqP, this is easily rewritten 

/oo pt 
dt P{t) / ds Vbe~^°(*~^)p(s) (64) 
-oo J —oo 



This is the main result of the paper: a class operator commuting with H, describing histories 
which enter the region A. It is valid in the approximation Eq.(62), which is sufficient to cover 
the key case of histories which, classically, intersect the boundary of A twice. 

For systems whose classical paths intersect the surfaces of interest more than two times, the 
semiclassical approximation Eq.(62) will not be valid, but the exact result Eq.(60) may still be 
used. It may also be of interest to explore higher order semiclassical approximations obtained 
by iterations of the basic result Eq.(58). This will not be explored here. 




Figure 4. A classical trajectory intersecting A enters at parameter time ta and leaves at 
parameter time t^. 



4. Properties of the Class Operators 

It is easy to show that the class operator for a single region Eq.(64) has a sensible classical limit. 
Classically, P{t) is a function on classical trajectories with P{t) = 1 when the classical trajectory 
is in A and is zero otherwise. Suppose a given trajectory q"'{t) enters A at some stage in its 
history so intersects the boundary twice (recalling we have essentially assumed no more than two 
intersections in the scmiclassical approximation Eq.(62)). For a given choice of parametrization 
of the trajectory, it enters A at parameter time ta and leaves at time ti, > ta (see Figure 4). (We 
may parameterize it in the opposite direction, with the same ultimate result). 
The derivative of P(s) is 

P{s) = 6is - ta) - 6{s - tb) (65) 
and we have ^ ^ 

Ca = [' dt f ds Foe-^°(*-') [5{s - ta) - 5{s - t^)] (66) 

*J ta ^ ta 

Since s < t < th, the second 5- function in Eq.(66) makes no contribution to the integral. This is 
exactly the desired property - the expression for Ca registers only the first intersection of the 
trajectory with the boundary, but not the second intersection. The integral is easily evaluated 
with the result, 

Ca = 1 - e-^°(*''-*°) (67) 
This is approximately 1, as required, as long as 

Voitb - ta) » 1 (68) 

We now give the broad picture in the quantum case and confirm some of the details in some 
simple models in the next section. Suppose we operate with Ca on an eigenstate \^x) of H. 
The time integrals in Eq.(64) are easily carried out and the result is 

CaI^a) = 27rVb<5(ii' - A) PGvP\^x) (69) 

where 

/•oo 

Gv = / e-^(^-^)*-^°* 
Jo 

—i 

H-X-iVo 



(70) 
(71) 



and we have used 

5{H-\) = — / dt e-^(-^-^)* (72) 

The P term is a current operator on the boundary E of A and may be written 

P = i[H,P] 

= -PnSi:{q) - 6j:{q)pn (73) 

where 

5^{q) = j^(P-\\q){q\ (74) 

is a 5-function operator in the surface S and Pn = rf'Pa is the component of the momentum 
operator normal to S with the outward pointing normal. The normal n" will depend on q in 
general so there may be an operator ordering issue in making p„ hermitian. The difference 
between different ordering will involve a term of the form Van", essentially the extrinsic 
curvature of S, and this will be small if A is sufficiently large and its boundary reasonably 
smooth. Note that it is sometimes also convenient to write P as 



1' 



<f'-Wq){Z-K){q\ (75) 



It is very useful to separate the current operator at the boundary into ingoing and outgoing 
parts according to the sign of Pn at the boundary 

P = {PU - {P)out (76) 
where (-P)m consists of the components of P with incoming momentum 

{P)in = -Pn0{-Pn)Si:iq) ' SEiq)Pn0{-Pn) (77) 

and similarly for {P)out- Examples of these expressions in particular models will be given in 
the following sections. The restriction to positive or negative p„ means that it is generally 
difficult to express (-P)i„ and {P)out in the form Eq.(75), involving the derivative pn = —idn, 
unless operating on a state (such as a WKB state) with simplifying properties. Note that these 
definitions require only that the local flux operator on a given surface can be split into ingoing 
and outgoing parts. To do so globally on a family of surfaces is generally impossible, essentially 
due to the problem of time, but fortunately this is not required here. 

The quantity Gy has the form of a Feynman propagator, Eq.(22), with Vq playing the role of 
the "ze prescription" , as long as Vq is sufficiently small (in comparison to an appropriate energy 
scale contained in H, such as po in the case of the Klein-Gordon equation). We therefore expect 
it to have properties similar to Gp, as described in Section ID. 

However, Vq is not set to zero exactly and this in fact means that Gy has two properties 
not possessed by Gp- First, the non-zero Vq produces a suppression for widely separated initial 
and final points. In a path integral representation of Gv{q,q'), the sum will be dominated 
by classical paths from q to q', to which one may associate the total parameter time r of the 
path. From the integral representation Eq.(70) it can be seen that Gy will have an overall 
exponential suppression factor of the form cxp( — Vqt). Thus, propagation with Gy will suppress 
configurations q and q' connected by classical trajectories of parameter time duration of greater 
than I/Vq. Second, recalling that the class operator is closely related to the 5- matrix for 
scattering off a complex potential there will be some reflection involved and this appears in 
properites of Gy - it will not exactly possess the Feynman properties of propagating positive 



frequencies to the future negative frequencies to the past, but may, for example, propagate some 
positive frequencies to the past. However, this "non-Feynman" propagation will be small for 
sufficiently small Vq, as we will see in the next section, so this possibility will be ignored. 
These properties of Gy are crucial to understanding the properties of Ca. First note that 

{q\PGvP\^x) = -iJ^ d^''(l' fA{q)Gviq, q'Wn - K)^xiq') (78) 

In this expression the propagator Gv{q,q') propagates from the boundary E to a point q on 
the interior. The fact that Gy is, approximately, a propagator of the Feynman type means two 
things. First, there will be initial and final points q, q' for which the propagator is very small 
(due to "propagation outside the lightcone", discussed in Section ID). Secondly, when it is not 
small, one would expect it to involve only ingoing modes at the boundary (to the extent that 
reflection is ignored). Hence the PGy terms effectively restrict the current operator P at the 
boundary to ingoing modes only - exactly the result we are looking for. We therefore replace 
P, with the current operator {P)in involving ingoing modes only 

Ca|*a) = 2'kVo5{H - A) PGy(P),„|*A) (79) 

Now note that since the PGy terms have done their job of selecting the ingoing modes, they 
may be eliminated. More precisely, we write P = 1 — P, where P is the projector onto the region 
outside A, and noting that that 

Vo5{H - \)Gv = S{H - A) (80) 

we find 

Ca\^x) = 27T5{H-X){P)in\^x) 

- 2Tr6{H-X)PGv{P)in\'i>x) (81) 

Consider the second term on the right-hand side. It consists of ingoing modes at the boundary 
of A propagated with Gy to a final point which is outside A. Semiclassically, this corresponds to 
paths which have to traverse the entire width of A, so if A is sufficiently large, this will take a long 
parameter time r and, as discussed above, Gy will contain a suppression factor of exp(— Vbr) 
in comparison to the first term in Eq.(81). This is exactly analogous to the suppression of the 
second term in the classical case, Eq.(67). We thus deduce that 

Ca|*a) ~ 27rS{H - X){P)in\^x) (82) 

Note that the approximation leading to dropping the second term in Eq.(81) also ensures that 
the result is independent of Vq, like the classical case in the appropriate limit. Eq.(82) is the 
main result of this section. 

Eq.(82) may also be expressed in terms of {P)out defined in Eq.(76). To see this, note that, 
since P = i[H, P], we have 

8{H - >!)P\^x) = i(A' - \)S{H - AOPI^-a) (83) 

This is zero for A = A', as long as (A[P[A) is well-defined, where |A) are eigenstates of H. (It 
may not be well-defined if P projects onto an infinite region - see Appendix A). Eq.(76) then 
implies 

5{H - A)(P)i„|*A) = 5{H - A)(P)„„t|*A) (84) 
so Eq.(82) may be expressed in terms of either the ingoing or outgoing flux or both. 



It is also useful to note that the right-hand side of Eq. (82) may be written 



2nSiH - \)iP)in\^x) = /s|*a) (85) 

where 

/oo 
dt e^^*(P)i„e-^^* (86) 
-oo 

This is an intersection number operator for ingoing flux at the boundary S of A. Hence the 
class operator is essentially the intersection number 1^,, and clearly commutes with H. It is 
classically equal to 1 for trajectories with ingoing flux at S (with no more than two intersections 
of the boundary, in the approximation we are using), and zero for trajectories not intersecting 
S. Of course, one may have guessed this approximate formula for the class operator, but a 
class operator is fundamentally defined as a product of projectors (or quasi-projectors) and 
the derivation given here makes it clear how this obvious guess arises from the fundamental 
definition. Note also the the class operator is hermitian in this case, so there is no need to 
consider the modified propagator Eq.(55). 

The form Eq.(82) may be used to check for decoherence of histories in specific models and 
we will see this later. When there is decoherence of histories, the probabilities are given by the 
average of a single class operator, Eq.(19), which in this case reads 

(*A'|Ca|^a) = 27r(^A|(i^)m|^A) '^(A - A') (87) 

Following the induced inner product prescription, we drop the (5-function on the right and then 
set A = A' = 0, so the probability in terms of a solution |^) of the Wheeler-DeWitt equation is 

{^\CAmphys = M^\iP)inm (88) 

This is essentially the ingoing Klein-Gordon flux on the boundary of A, as expected. (The factor 
of 27r relates to the induced inner product as described in Appendix A of Ref.[16]). 

The above derivation is verified in more detail in Ref.[16] using specific examples. 
Furthermore, it is also shown there that the form of the class operator for n regions is, 

Cn = —. (-fsi-^Ea ■ ■ ■ -^E„ + permutations) (89) 
n! 

where "permutations" means add all possible permutations of the n regions, to give a total of 
n! terms. This final result has a particularly natural form which also suggests that it may have 
simple path integral representations. This will be explored elsewhere. (See also Ref.[42] for 
similar expressions, derived using a detector model for quantum cosmology). 

5. WKB Regime 

The most important case in which to check the ideas developed above is in the WKB regime. 
In the oscillatory regime, the solutions to the Wheeler-DeWitt equation have the form 

= Re'^ (90) 

where R and S obey Eqs.(24), (25) as described earlier. More generally, the wave function is 
a superposition of WKB wave functions but we consider first the case of a single term. The 
heuristic interpretation of such states was described in Section ID. Our aim is to show that 
the decoherent histories analysis reproduces the heuristic scheme. WKB states are locally plane 
wave states of the type considered in the previous Section, so we will appeal to that analysis to 
understand the properties WKB states. 



We consider the action of the class operator for a single region A on a WKB state. We first 
consider the class operator Eq.(82) acting on WKB state (regularized by making it an eigenstate 
of H with eigenvalue A) 

(9|Ca|*a) = 2'K{q\5{H-\){P)in\^x) 

= 2'Ki j cr-\' {q\5{H - \)\q'){t - K)^xiq') (91) 

Here Sj„ denotes the sections of the boundary where the flux is ingoing, 

n • VS" < (92) 

where is the outward pointing normal. The key property of the WKB wave function ^x{q') 
is that it has momentum p = VS'(g') at each point q' on the boundary. When the operator 
6{H — A) is applied, from the representation Eq.(72), we see that it's effect is to evolve the 
state '^x{q') (restricted to Sj„) forwards and backwards in parameter time, and then integrates 
over all times. Semiclassically, the evolution of the state "^x{<l') will be concentrated along the 
classical trajectories defined by initial positions q' in Sjjj and momenta p — ViS(5'), with some 
spreading of the wave packet, but this will be small if A is reasonably large. (This is analogous 
to the model of Section 5B). 

We thus see the following: the wave function (q'|Ca|^'A)) in Eq.(91) is spatially localized 
around the tube of classical trajectories passing through A with momenta p = VS (depicted in 
Figure 2). This may be approximately written in the alternative form 

{q''\CA\^x)^0{TA'e) Re'^ (93) 

where e > is a small parameter to regularize the ^-function at zero argument. Here, T^iq) is 
the parameter time spent by the classical trajectory qci{t) (with initial value q and momentum 
p = VS{q)) in the region A and may be written 

/oo 
dt fA{qcl{t)) (94) 
-oo 

This has the property that 

V5 • Vta = (95) 

since VS-V simply translates along the classical trajectories. It follows from Eq.(25) that Eq.(93) 
is in fact a WKB solution to the Wheeler-DeWitt equation, since the ^-function essentially 
modifies the prefactor R but in a way that it still satisfies Eq.(25). (Of course, this is expected 
because Ca commutes with H). 

Eq. (93) is a very useful result and allows us to check for decoherence very easily. The action 
of the class operator for not entering A is clearly 

(cz1CA|^A)~^(e-TA) i?e^^ (96) 

which is a WKB state localized on the set of trajectories not entering A. It immediately follows 
that 

(^a'IC'aCaI^'a) ~0 (97) 

since the two states Eqs.(93), (96) are localized about complementary regions. There is therefore 
approximate decoherence of histories for a single WKB packet and for histories entering or not 
entering a single region A, as long as A is sufficiently large. 



The key reason for the decoherence with a single WKB packet is related to the approximate 
determinism of the WKB wave functions: fixing values of position to lie on Xj^n also fixes the 
momenta, since p = VS", so that Eq.(91) is concentrated along a tube of classical trajectories. 

More generally, the initial state will be a superposition of WKB wave packets, 

^ = Y^ Rke'^^ (98) 
k 

The component states in this sum are typically approximat ely orthogonal to each other as long 
as the phases S'fc are sufficiently different. (This will depend on the detailed dynamics of the 
model). Following the analogous example in the simple models of Section 5B, we would expect 
that the class operators will not disturb the approximate orthogonality of these states as long 
as the region A is sufficiently large - we expect the cross terms in the decoherence functional to 
average to zero because of oscillations. Therefore, superpositions may in practice be treated as 
mixtures at sufficiently coarse-grained scales. Note that this statement also applies the special 
state, 

m = R (e*^ + e"^^) (99) 

which arises from the Hartle-Hawking "no boundary" proposal [45]. That is, the interference 
between the two terms may be neglected. 

Given decoherence, the probabilities are given by the general expression Eq.(88). It then 
follows that the probabilities for entering A coincide with Eq.(29), the sought-after result. 

Now consider the case of probabilities for histories entering two regions, as described by the 
class operator Eq.(89). Since the two-region class operator is a sum of products of one-region 
class operators, its effect on the WKB wave functions is easy to see. The action of a single class 
operator gives Eq.(93). But since this is still a wave function of the WKB type, the action of a 
second class operator yields 

('zIC'aiaJ^a) ~ ^(tai - e)^(rA2 - e) Re'^ (100) 

That is, it is a WKB wave function but restricted in such as way that it's flux passes through 
both regions. See Figure 5. It is again easy to see that there is decoherence of histories and the 
probability is given by an expression of the form Eq.(29), but with the integral over the subset 
of incoming flux at Ai which goes on to intersect A2. 

From the above observations one can also see why the exponential WKB wave functions will 
not lead to decoherence of histories. The exponential wave functions have the form 

^ = Re-^ (101) 

where / is real. The key difference between states of this type and the oscillatory type 
Eq.(90) is that they do not have a correlation between positions and momenta [29]. One 
would therefore expect the evolution of the state Eq.(91) to be spread all over the configuration 
space and not concentrated around a particular region. That is, these states do not have the 
approximate determinism of the oscillatory states. The states Ca|^) and CaI^) would then 
not be approximately orthogonal so there will be no decoherence of histories. 

The decoherence of histories described here has arisen because of the approximate 
determinism of the oscillatory WKB states, together with the approximate orthogonality 
properties that arise when the regions A are sufficiently large. At finer grained scales, 
decoherence of histories may only be possible in more complicated models in which there is 
an environment of some sort. Models along these lines, in more basic approaches to quantum 
cosmology, have been considered previously [46], and one might expect that they may be adapted 
to the decoherent histories approach to quantum cosmology. (See also Ref.[ll]). 

In summary, we have derived from the decoherent histories approach the probabilities 
normally used in the heuristic WKB interpretation. To address these issues in more detail 
will require more specific quantum cosmological models. This will be considered elsewhere. 
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Figure 5. The class operator CaiA2 for two regions operating on a WKB wave function 
produces another WKB state locaUzed around the flux passing through both regions. This is 
the same as the flux entering Ai across the surface S12. 

6. Discussion and Further Issues 

We have presented a properly defined quantization procedure for quantum cosmology using the 
decoherent histories approach to quantum theory and derived from this the frequently used but 
heuristic WKB interpretation, involving fluxes of the WKB wave function. 

The key idea was to use a complex potential to define the class operators for not entering a 
region of configuration space. This method is adequately justified by its successful use in the 
arrival time problem in non-relativistic quantum theory. We showed that the class operators 
defined in this way have all the desired properties - they have the correct classical limit, are 
compatible with the constraint equation, and do not have difficulties with the Zeno effect. In 
a semiclassical approximation, they have an appealing form in terms of intersection number 
operators. They give sensible results in simple models and there is approximate decoherence of 
histories for certain types of initial state at sufficiently coarse grained scales. 

Future papers will address the more detailed application of this approach to specific models 
and will also undertake a comparison of the decoherent histories approach described here to 
other approaches [5]. 
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